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DATA  AND  NOTATION 


Model  Data:  (See  Fig.  1) 

1.  The  model  flies  3  kg  of  its  own  weight  at  35  m/sec. 

2.  Overall  model  lift  =  0.875  kg/taileron  degree  at  35  m/sec. 

3.  Span  =  1.8  m. 

4.  Length  =  2.0  m 

Nooenclature: 

Kact  Taileron  actuator  static  gain  constant. 

Kj  Pitch  rate  weighting. 

K2  Pitch  attitude  weighting. 

K3  CG  displacement  weighting. 

K4  CG  velocity  weighting. 

ly  Pitch  inertia  about  pivot 

(Initially  23.3  kg.m^,  modified  to  21.2  kg.m^). 


Aircraft  mass 

(Initially  66.6  kg.  modified  to  62.8  kg) 

Reference  area  (0.6  m^). 

Laplace  transform  variable. 

Vertical  velocity  at  CG,  positive  downwards  (m/sec). 

Tunnel  velocity  (m/sec). 

(Initially  1.441m) 

Distance  to  GG  (Modified  to  1.543m)  Positive  to  rear, 

measured  from  an 

Distance  to  pivot  (Initially  1.51m)  origin  near  the  nose 

(Modified  to  1.55m) 

Distance  from  mass  relief  spring  to  CG,  positive  for  the  CG  in 
front  of  the  spring. 

(Initially  0.068m,  modified  to  0.007m). 

Pitch  rate,  positive  nose  up  (rad/sec). 

Taileron  angle  (rad). 

Taileron  trim  angle  (measured  at  -5  deg  at  35  m/sec). 

Total  taileron  angle  available  (i  6  degrees). 
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Pitch  angle,  positive  nose  up  (rads). 
Density  of  air  (1.225  kg/m^). 


1^ 

mi 


.s  .'v'  * 
rX-Vv: 


V  Vv* 


Taileron  actuator  time  constant  (0.01  sec). 

Time  constant  of  pseudo-integrator  (6.0  sec). 

Time  constant  of  sensor  noise  filter.  (0.03  sec). 

Natural  frequency  of  pitch  spring  (1.885  rad/sec). 

Natural  frequency  of  heave  spring  (4.712  rad/sec). 

Assumed  damping  of  pitch  spring  (2%). 

Assumed  damping  of  heave  spring  (5%). 

Vertical  displacement  at  CG,  positive  downwards  (m). 

Wing  lift  coefficient  (-  2.902/rad,  acting  at  x  =  1.657m). 

Tail  lift  coefficient  (-  0.672/rad,  acting  at  x  =  2.028m). 

Coefficient  of  lift  due  to  pitch  rate  (-  0.581/rad  acting  at 
X  =  2.072m) 


NOTE There  is  a  hidden  (unit)  length  in  the  Zq  term  due  to 
Kaynes'  method  of  derivation  of  the  coefficients  (Ref.  1). 

Denotes  estimated  value. 

Dot  over  a  symbol  denotes  a  derivative  with  respect  to  time. 


1.  INTRODUCTION 


The  author  was  attached  to  Flight  Systems  Department,  RAE, 
Farnborough  UK,  from  December  1982  to  February  1984.  Part  of  that  time 
was  spent  on  a  task  for  Structures  and  Materials  Department,  under  the 
direction  of  Mr.  Colin  Skingle  of  the  Experimental  Dynamics  Group.  This 
group  was  involved  in  a  collaborative  programme,  under  the  auspices  of 
GARTEur,  on  gust  load  alleviation  and  flutter  suppression  research.  The 
programme  was  divided  into  three  phases,  classified  as  follows  (Ref.  2);- 

Phase  1:  Production  and  demonstration  of  flutter  suppression 

control  laws  on  a  wind  tunnel  aircraft  model  up  to  1.7 
Qp,  where  Qp  is  the  flutter  dynamic  pressure. 

Phase  2:  Demonstration  of  flutter  suppression  control  laws  to 

1.5  Qp  in  the  presence  of  gust  inputs  -  both  Dryden 
spectra  and  discrete  gusts  being  programmable  from  a 
gust  generator  upstream  of  the  model.  My  activities 
were  confined  to  this  phase. 

Phase  3:  Demonstration  of  flutter  suppression  and  gust  load 

alleviation  (GLA)  laws  simultaneously,  to  1.5  Qp, 
using  the  gust  generator. 


The  participating  countries  were  UK  (RAE/BAe),  Holland  (NLR),  W. 
Germany  (DFVLR/MBB)  and  France,  using  the  wind  tunnel  facilities  at  DFVLR, 
Gottingen  ,  W.  Germany  and  a  model  supplied  by  DFVLR/MBB.  Phase  1  was 
completed  in  1982,  Phase  2  completed  in  November  1983  and  Phase  3  was 
scheduled  for  1984. 

2.  SIMULATION  STUDIES 

The  modi fi ed-Tornado  model  was  mounted  so  that  it  was  constrained  to 

heave  and  pitch  freedoms  only  (Fig.  1).  The  high  pressure  hydraulic 

supply  to  the  model  for  control  surface  actuation  was  through  looped 

flexible  metal  tubes,  arranged  like  a  safety  pin,  with  low  torsional  and 
vertical  stiffness.  Phase  1  trials  had  shown  that  the  rigid-body  heave 
and  pitch  modes  were  very  lightly  damped.  One  of  my  activities  in 
preparation  for  the  Phase  2  trials  was  to  investigate  the  rigid-body 

response  of  the  model  and  to  develop  a  stability  augmentation  system  (SAS) 
control  law  for  it.  The  RAE  designed  control  system  simulation  and 

analysis  computer  package  'TSIM'  (Refs.  3,  4)  was  used  to  program  the 

rigid  body  and  control  law  model  simulation  (Appendices  A,  B). 

The  roots  of  the  author's  rigid-body  model  modes  were  compared  with 
those  of  a  complete  rigid-body  plus  aeroelastic  model  developed  by  Ian 
Kaynes  of  Structures  Dept.  RAE  in  order  to  verify  the  simpler  model.  The 

comparisons  were  performed  for  two  values  of  with  a  variety  of  tunnel 

speeds.  These  comparisons  are  depicted  in  Tables  A  and  B  and  show  that 
results  of  the  simpler  model.  Type  2,  are  not  significantly  different 
from  those  of  the  full  model.  Type  1. 


The  root  locus  of  the  basic  aircraft  with  changing  velocity  V  is 
depicted  in  Fig.  2.  for  nominal  values  of  =  0.068m  and  Xg  =  1.441m. 

The  unstable  rigid-body  mode  is  apparent  at  V  =  28  m/sec  although  in 
practice,  the  damping  was  such  that  this  mode  was  marginally  stable  over 
the  whole  speed  range. 

Phase  1  trials  also  showed  that  with  the  model  trimmed  at  nominal  Xp 
and  Xg  values  and  Xgpp  =  0.007m  (i.e.  mass  relief  spring  approximately  at 
the  pivot),  the  tailplane  trim  angle  was  found  to  be  large  (approximately 
-5  degrees).  Most  of  the  available  ±6°  movement  was  used  just  for 
trimming  and  very  little  was  left  for  using  the  tail  as  an  active 
control.  Suggested  remedies  were;- 

(a)  to  move  the  CG  rearward  and  hence  reduce  the  required  tail 
download  to  trim,  and/or 

(b)  to  move  the  spring  attachment  forward. 

Both  these  options  were  explored  using  TSIM  and  both  tended  to 
destabilize  the  rigid  aircraft  modes  even  further,  as  evident  in  the  root 
locus  diagrams.  Figs  3  and  4. 

Therefore,  the  stability  augmentation  system  (SAS)  was  required  to 
provide  basic  stability  and  if  possible,  to  extend  the  damping  of  these 
modes  to  a  value  more  representative  of  that  of  free  flying  aircraft. 
Relevant  feedback  sensors  available  were  accelerometers  in  the  fuselage 
(nose  and  near  the  CG)  and  a  pitch  rate  gyro. 
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As  a  design  exercise,  it  was  assumed  that  the  designer  had  latitude 
in  the  values  of  X^pfj  and  Xg.  Thus  one  design  goal  was  to  choose  values 
of  Xgpp  and  Xg  for  optimum  'basic  aircraft'  modal  damping  before  applying 
the  SAS  control  law.  However  as  the  study  progressed,  information  was 
received  from  DFVLR,  Gottingen  that  the  spring  and  CG  positions  were 
immutable.  The  simulation  studies  showed  that  when  the  SAS  law  was 
applied,  root  sensitivity  to  variations  in  spring  attachment  position  and 
CG  position  became  small.  Therefore  even  if  these  positions  had  been 
variable,  there  would  have  been  little  value  in  changing  them. 

The  root  locus  of  the  basic  aircraft  with  changing  velocity  V  is 
depicted  in  Fig.  2  for  the  original  values  of  X^pp  =  0.068m  and  Xg  = 

1.441m.  The  unstable  rigid  body  mode  is  apparent  at  V  =  28  m/sec.  Fig.  3 

shows  the  root  locus  for  a  fixed  speed  V  =  35.0  m/sec  with  varying  X^pp. 
The  stabilizing  effect  of  increasing  X^pp  (i.e.  moving  the  location  of  the 
spring  forward)  is  clearly  shown.  The  similar  stabilizing  effect  of 
moving  the  CG  location  further  aft  (increasing  Xg)  is  shown  in  Fig.  4. 
These  effects  are  summarised  in  the  dual -parameter  families  of  root  loci 
in  Figs.  5  and  6.  From  these  plots,  values  of  Xgpp  =  0.068m  and  Xg  =  1.6m 
were  chosen  to  give  optimum  open-loop  stability  of  the  basic  aircraft 
without  feedback.  This  reference  open-loop  condition  (Fig.  7)  was  used 
for  initial  studies  of  an  SAS  control  law.  Fig.  7  shows  that  the  damping 
ratio  of  the  dominant  open-loop  poles  for  the  values  chosen  is  0.15  at  35 

m/sec.  The  root  sensitivity  to  Xgpp  is  quite  large  in  this  region,  as 

shown  in  the  time  responses  to  a  taileron  step  for  Xgpp  values  of  0.07, 
0.08  and  0.09m  in  Figs.  8,  9,  and  10.  One  root  has  become  positive  real 
in  Fig.  10  for  a  0.01m  increase  in  X5pp. 
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3. 


CONTROL  LAW  DESIGN 


The  four  states  corresponding  to  the  tv/o  Rigid  Body  Modes  (RBM)  of 
pitch  and  heave,  i.e.  pitch  rate  q,  pitch  angle  6  ,  CG  vertical 
displacement  Zq  and  velocity  V2g  were  fed  back  according  to  the  law 

S  =  (K^  q  +  K2  0  +  Zg  +  V^q).  G^ 

^ACT 

where  G.  - - ,  the  actuator  transfer  function. 

»  *  ^ACT-=) 

Numerous  root-loci  and  time  responses  were  evaluated  using  trial 
weightings  of  the  four  states  before  the  following  RBM  state  feedback  law 
was  obtained 


6  =  (0.4  q  +  9  +  0.2  Zq).G^  for  Xg  =  1.6m 

and  ^  0.07m 

The  root  locus  of  this  law  for  varying  wind  speed  is  shown  in  Fig. 
11(a)  using  a  first-order  model  for  the  taileron  actuator,  with  time 
constant  of  0.01  seconds. 


(6) 


For  practical  implementation  of  this  law,  0  and  Zq  were  estimated  as 


fol lows 


since 


6  =  q 


q 

6  =  - 
s 


Using  pseudo-integration,  the  integral  time  constant  Tjfj-j-  was 
selected  to  be  about  10  times  the  time  constant  of  the  mode  to  be 
controlled.  The  pseudo-integrator  was  also  gain-matched  (for  large  s) 
with  the  perfect  integrator. 
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This  yielded  0  =  * - r 

t  TNT^ 


Similar  reasoning  on  the  estimates  and  Zg  yielded 


(7) 


Thus  the  block  diagrain  of  Fig.  12  is  obtained,  with  the  addition  of 
a  sensor  noise  filter  and  first-order  model  of  the  actuator,  with  the 
estimate  6  given  by 
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The  root  locus  of  the  estimated  RBM  state  feedback  is  shown  in  Fig. 
11(b)  for  the  same  weightings  determined  previously,  i.e.  Ki  =  0.4,  K2  = 
1.0,  K3  =  0.2  and  =  0.  The  comparison  with  Fig.  11(a)  shows  that  the 
behaviour  of  the  pitch-and  heave-mode  roots  is  essentially  similar  with 
estimated  state  feedback. 


The  TSIM  models  of  the  'basic  aircraft',  aircraft  with  RBM  state 
feedback,  and  aircraft  with  estimated  RBM  state  feedback  are  given  in 
Appendix  B. 


Subsequently,  information  from  Gottingen  on  the  final  fixed 
positions  of  and  Xg  were  obtained  (0.007m  and  1.543m  respectively) 

together  with  updated  values  of  M  and  1^. 
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The  basic  aircraft  root  locus  for  the  updated  data  is  shown  in  Fig. 
13,  where  the  aircraft  becomes  unstable  at  24  m/sec  and  stays  unstable 
through  to  60  m/sec.  (Compare  with  Fig.  2  where  stability  is 
re-established  at  about  55  m/sec.) 

A  similar  study  to  that  previously  undertaken  was  performed  to 
choose  an  acceptable  control  law.  Because  this  was  to  be  the  final 
version  for  practical  implementation,  a  simplified  pitch-rate  feedback  law 
was  evaluated  and  subsequently  judged  acceptable.  The  root  locus  of  the 
pitch-rate  feedback  law  is  shown  in  Fig.  14.  The  time  responses  for  a 
taileron  step  displacement  of  5  degrees  with  three  wind  speeds  are  shown 
in  Fig.  15.  A  similar  time  response  for  the  best  evaluated  RBM  state 
feedback  law,  in  this  case  given  by  - 


6  =  0.5  q  +  Q  +  0*1  Zq  -  0-1 

is  shown  in  Fig.  16.  This  shows  a  reduction  in  an  already  small  heave 
amplitude  of  about  10mm  at  the  expense  of  slightly  degraded  pitch 
damping.  Thus  it  is  apparent  there  is  little  benefit  from  feeding  back 
all  the  rigid  body  modes  in  this  case.  In  practice,  the  law  ^  =  0.5  q  was 
used  with  success. 

The  time  response  at  35  m/sec  for  the  modified  data  'basic  aircraft' 
is  also  shown  in  Fig.  17  where  ±  1  metre  heave  amplitudes  are  seen  to 
build  up  within  5  seconds. 


This  divergent  rigid  body  oscillation  was  the  mechanism  which 
slammed  the  model  into  the  end  stops  during  the  French  team's  trials. 
Fortunately,  the  damage  was  repaired  quickly  and  did  not  significantly 
affect  the  performance  of  the  subsequent  RAE  team. 

4.  CONCLUSION 

This  document  has  outlined  some  work  undertaken  for  Structures  and 
Materials  Department,  RAE  as  a  small  part  of  the  author's  attachment. 
Because  the  task  was  also  a  valuable  learning  exercise,  the  study  of  a 
suitable  feedback  law  was  made  more  comprehensively  than  normally 
justifiable  on  a  simple  model. 
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APPENDIX  A  (CONT'O) 


Similarly  the  Moment  eqn.  is  - 


ly.q  =  PSV 


(1.657  -  Xg)( 


Vy  +e)  +  PSVZq. 


-  .  (Zg  +  Xgpj^  .0).  Xgpp 


-  (Zg  +  Xjpp.q)  .  Xgpp 


-  0  -  2  ly  Ce  «e  .  q 

+  PSV^Z,,,  (2.028  -  X.)  .  <5 


(2.0723 


TSIM  Nomenclature 
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APPENDIX  B  {CONT*D) 

B2,  TSIN  ■odel  of  basic  aircraft  rigid-body  response 
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APPENDIX  B  (CONT'D) 

B3.  Aircraft  rigid  body  wodcl  with  full  state  feedbacfc  \am 
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DERIVATIVE  ZWT , XG , lY , FREQTH , DAMPTH 
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